Stability and energy budget of pressure-driven collapsible channel flows

Although self-excited oscillations in collapsible channel flows have been extensively studied, our understanding of their origins and mechanisms is still far from complete. In the present paper, we focus on the stability and energy budget of collapsible channel flows using a fluid–beam model with the pressure-driven (inlet pressure specified) condition, and highlight its differences to the flow-driven (i.e. inlet flow specified) system. The numerical finite element scheme used is a spine-based arbitrary Lagrangian–Eulerian method, which is shown to satisfy the geometric conservation law exactly. We find that the stability structure for the pressure-driven system is not a cascade as in the flow-driven case, and the mode-2 instability is no longer the primary onset of the self-excited oscillations. Instead, mode-1 instability becomes the dominating unstable mode. The mode-2 neutral curve is found to be completely enclosed by the mode-1 neutral curve in the pressure drop and wall stiffness space; hence no purely mode-2 unstable solutions exist in the parameter space investigated. By analysing the energy budgets at the neutrally stable points, we can confirm that in the high-wall-tension region (on the upper branch of the mode-1 neutral curve), the stability mechanism is the same as proposed by Jensen and Heil. Namely, self-excited oscillations can grow by extracting kinetic energy from the mean flow, with exactly two-thirds of the net kinetic energy flux dissipated by the oscillations and the remainder balanced by increased dissipation in the mean flow. However, this mechanism cannot explain the energy budget for solutions along the lower branch of the mode-1 neutral curve where greater wall deformation occurs. Nor can it explain the energy budget for the mode-2 neutral oscillations, where the unsteady pressure drop is strongly influenced by the severely collapsed wall, with stronger Bernoulli effects and flow separations. It is clear that more work is required to understand the physical mechanisms operating in different regions of the parameter space, and for different boundary conditions

Study of cardiovascular function using a coupled left ventricle and systemic circulation model

To gain insight into cardio-arterial interactions, a coupled left ventricle-systemic artery (LV–SA) model is developed that incorporates a three-dimensional finite-strain left ventricle (LV), and a physiologically-based one-dimensional model for the systemic arteries (SA). The coupling of the LV model and the SA model is achieved by matching the pressure and the flow rate at the aortic root, i.e. the SA model feeds back the pressure as a boundary condition to the LV model, and the aortic flow rate from the LV model is used as the input for the SA model. The governing equations of the coupled system are solved using a combined immersed-boundary finite-element (IB/FE) method and a Lax–Wendroff scheme. A baseline case using physiological measurements of healthy subjects, and four exemplar cases based on different physiological and pathological scenarios are studied using the LV–SA model. The results of the baseline case agree well with published experimental data. The four exemplar cases predict varied pathological responses of the cardiovascular system, which are also supported by clinical observations. The new model can be used to gain insight into cardio-arterial interactions across a range of clinical applications

Constricted channel flow with different cross-section shapes

Pressure driven steady flow through a uniform circular channel containing a constricted portion is a common problem considering physiological flows such as underlying human speech sound production. The influence of the constriction’s cross-section shape (circle, ellipse, circular sector) on the flow within and downstream from the constriction is experimentally quantified. An analytical boundary layer flow model is proposed which takes into account the hydraulic diameter of the cross-section shape. Comparison of the model outcome with experimental and three-dimensional numerically simulated flow data shows that the pressure distribution within the constriction can be modeled accurately so that the model is of interest for analytical models of fluid–structure interaction without the assumption of two-dimensional flow

A modified Holzapfel-Ogden law for a residually stressed finite strain model of the human left ventricle in diastole

In this work, we introduce a modified Holzapfel-Ogden hyperelastic constitutive model for ventricular myocardium that accounts for residual stresses, and we investigate the effects of residual stresses in diastole using a magnetic resonance imaging–derived model of the human left ventricle (LV). We adopt an invariant-based constitutive modelling approach and treat the left ventricular myocardium as a non-homogeneous, fibre-reinforced, incompressible material. Because in vivo images provide the configuration of the LV in a loaded state even in diastole, an inverse analysis is used to determine the corresponding unloaded reference configuration. The residual stress in this unloaded state is estimated by two different methods. One is based on three-dimensional strain measurements in a local region of the canine LV, and the other uses the opening angle method for a cylindrical tube. We find that including residual stress in the model changes the stress distributions across the myocardium and that whereas both methods yield qualitatively similar changes, there are quantitative differences between the two approaches. Although the effects of residual stresses are relatively small in diastole, the model can be extended to explore the full impact of residual stress on LV mechanical behaviour for the whole cardiac cycle as more experimental data become available. In addition, although not considered here, residual stresses may also play a larger role in models that account for tissue growth and remodelling

Correlation of mechanical factors and gallbladder pain

Acalculous biliary pain occurs in patients with no gallstones, but is similar to that experienced by patients with gallstones. Surgical removal of the gallbladder (GB) in these patients is only successful in providing relief of symptoms to about half of those operated on, so a reliable pain-prediction model is needed. In this paper, a mechanical model is developed for the human biliary system during the emptying phase, based on a clinical test in which GB volume changes are measured in response to a standard stimulus and a recorded pain profile. The model can describe the bile emptying behaviour, the flow resistance in the biliary ducts, the peak total stress, including the passive and active stresses experienced by the GB during emptying. This model is used to explore the potential link between GB pain and mechanical factors. It is found that the peak total normal stress may be used as an effective pain indicator for GB pain. When this model is applied to clinical data of volume changes due to Cholecystokinin stimulation and pain from 37 patients, it shows a promising success rate of 88.2% in positive pain prediction

Volumetric growth of soft tissues evaluated in the current configuration

The growth and remodelling of soft tissues plays a significant role in many physiological applications, particularly in understanding and managing many diseases. A commonly used approach for soft tissue growth and remodelling is volumetric growth theory, introduced in the framework of finite elasticity. In such an approach, the total deformation gradient tensor is decomposed so that the elastic and growth tensors can be studied separately. A critical element in this approach is to determine the growth tensor and its evolution with time. Most existing volumetric growth theories define the growth tensor in the reference (natural) configuration, which does not reflect the continuous adaptation processes of soft tissues under the current configuration. In a few studies where growth from a loaded configuration was considered, simplifying assumptions, such as compatible deformation or geometric symmetries, were introduced. In this work, we propose a new volumetric growth law that depends on fields evaluated in the current configuration, which is residually stressed and loaded, without any geometrical restrictions. We illustrate our idea using a simplified left ventricle model, which admits inhomogeneous growth in the current configuration. We compare the residual stress distribution of our approach with the traditional volumetric growth theory, that assumes growth occurring from the natural reference configuration. We show that the proposed framework leads to qualitative agreements with experimental measurements. Furthermore, using a cylindrical model, we find an incompatibility index that explains the differences between the two approaches in more depth. We also demonstrate that results from both approaches reach the same steady solution published previously at the limit of a saturated growth. Although we used a left ventricle model as an example, our theory is applicable in modelling the volumetric growth of general soft tissues

On the chordae structure and dynamic behaviour of the mitral valve

We develop a fluid-structure interaction (FSI) model of the mitral valve (MV) that uses an anatomically and physiologically realistic description of the MV leaflets and chordae tendineae. Three different chordae models — complex, “pseudo-fibre”, and simplified chordae — are compared to determine how different chordae representations affect the dynamics of the MV. The leaflets and chordae are modelled as fibre-reinforced hyperelastic materials, and FSI is modelled using an immersed boundary-finite element (IB/FE) method. The MV model is first verified under static boundary conditions against the commercial FE software ABAQUS, and then used to simulate MV dynamics under physiological pressure conditions. Interesting flow patterns and vortex formulation are observed in all three cases. To quantify the highly complex system behaviour resulting from FSI, an energy budget analysis of the coupled MV FSI model is performed. Results show that the complex and pseudo-fibre chordae models yield good valve closure during systole, but that the simplified chordae model leads to poorer leaflet coaptation and an unrealistic bulge in the anterior leaflet belly. An energy budget analysis shows that the MV models with complex and pseudo-fibre chordae have similar energy distribution patterns, but the MV model with the simplified chordae consumes more energy, especially during valve closing and opening. We find that the complex chordae and pseudo-fibre chordae have similar impact on the overall MV function, but that the simplified chordae representation is less accurate. Because a pseudo-fibre chordal structure is easier to construct and less computationally intensive, it may be a good candidate for modelling MV dynamics or interaction between the MV and heart in patient-specific applications

Preparation and quality evaluation of Ibuprofen and Piphenhydramine hydrochloride orally disintegrating tablets

The present study was designed to evaluate the preparation and quality of ibuprofen and diphenhydramine hydrochloride orally disintegrating tablets. The formulation was optimized with disintegration time and taste as reference parameters by single-factor and orthogonal experiments. Then the hardness, disintegrating time and dissolution were examined. The best prescription was found to contain: ibuprofen 200 mg, diphenhydramine HCl 25 mg, citric acid 30 mg, microcrystalline cellulose 115 mg, mannitol 48.5 mg, polyvinylpolypyrrolidone 40 mg, sodium lauryl sulphate 10 mg, magnesium Stearate 3.5 mg, gum arabic 1.5 mg, aspartame 15 mg, steviosin 1.5 mg, and sodium bicarbonate 10 mg. The orally disintegrating tablets were congruent and smooth with pleasant taste in mouth. The disintegration time was less than 60 s, and more than 80 % of ibuprofen and diphenhydramine HCl were dissolved within 10 min. Consequently, the formulation design is reasonable and the process of preparation is feasible.Colegio de Farmacéuticos de la Provincia de Buenos Aire

Anisotropic behaviour of human gallbladder walls

Inverse estimation of biomechanical parameters of soft tissues from non-invasive measurements has clinical significance in patient-specific modelling and disease diagnosis. In this paper, we propose a fully nonlinear approach to estimate the mechanical properties of the human gallbladder wall muscles from in vivo ultrasound images. The iteration method consists of a forward approach, in which the constitutive equation is based on a modified Hozapfel–Gasser–Ogden law initially developed for arteries. Five constitutive parameters describing the two orthogonal families of fibres and the matrix material are determined by comparing the computed displacements with medical images. The optimisation process is carried out using the MATLAB toolbox, a Python code, and the ABAQUS solver. The proposed method is validated with published artery data and subsequently applied to ten human gallbladder samples. Results show that the human gallbladder wall is anisotropic during the passive refilling phase, and that the peak stress is 1.6 times greater than that calculated using linear mechanics. This discrepancy arises because the wall thickness reduces by 1.6 times during the deformation, which is not predicted by conventional linear elasticity. If the change of wall thickness is accounted for, then the linear model can used to predict the gallbladder stress and its correlation with pain. This work provides further understanding of the nonlinear characteristics of human gallbladder